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ABSTRACT 

We have cross-correlated the SDSS DR3 [Schneider et all (|2005D quasar catalog with the XMM- 
Newton archive. Color and redshift selections {g — r > 0.5 and 0.9 < z < 2.1) result in a sample of 
17 red, moderate redshift quasars. The redshift selection minimizes possible contamination due to 
host galaxy emission and Lya forest absorption. Both optical and X-ray information are required to 
distinguish between the two likely remaining causes of the red colors: 1) dust-reddening and 2) an 
intrinsically red continuum. We find that 7 of 17 quasars can be classified as probable 'intrinsically 
red' objects. These 7 quasars have unusually broad Mgll emission lines (<FWHM> = 10,500 km 
s^^), moderately flat, but unabsorbed X-ray spectra (< T >=1.66±0.08), and low accretion rates 
{M / Msdd ^ 0.01). We suggest low accretion rates as a a possible physical explanation for quasars 
with intrinsically red optical continua. We find that 8 of 17 quasars can be classified as dust-reddened. 
Three of these have upper-limits on the absorption column from X-ray spectral fits of Nh = 3-13 x 10^^ 
cm^, while the other five quasars must be absorbed by at least Nh — cm^ in order to be consistent 
with a comparably selected aox — luv distribution. Two objects in the sample are unclassified. 
Subject headings: accretion disks — galaxies: active — quasars: general 



1. INTRODUCTION 

Until recently, quasars have typically been selected as 
point-like objects with colors bluer than stars (e.g. U — B 
< -0.44, Schmidt & Green 1983). This limited the possi- 
ble spectral energy diagrams (SEDs) that known quasars 
could have. A composite SED shows that quasars' blue 
colors are due to the Big Blue Bump (BBB) continuum 
feature that dominates the optical and U V spectrum 
(jMalkan fc Sargent 1 119821 : lElvis et all Il994[ ) . However, 
surveys using different selection mechanisms discovered 
that quasars can have red colors as well, indicating a 
diminished or missing BBB. While in many cases, this 
reddening is due to dust obscuration, some cases may be 
due to changes in the accretion disk thought to power the 
optical/UV continuum. We identify such a population in 
this paper. 
Red quasars were first discovered in the 



1995 ) sample of radio-loud quasars. 
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19951 ) suggested dust-reddening as the cause of the red 
Bj — K colors and estimated that red colors could cause 
as much as 80% of the quasar population to go unde- 
tected in existing optical surveys. Some 200 members of 
a new population of red, radio-quiet quasars were discov- 
ered in the near-infrared (1-2 ^,m) with th e Two Micron 
All Sky Survey (2MASS) by iCutri et al.l ((200l . The 
space density of the red 2MASS quasars was found to be 
at least equ a l to th at of optical and UV-selected quasars. 
ICutri et all (|200^ also suggested dust-reddening to ex- 
plain the red colors of the 2MASS red quasars because 
the near-IR color distribution is consistent with a red- 
dening of Ay = 1-5 magnitudes. In addition, ^^10% 
of the predominantly radio-quiet sample is highly po- 
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larized (P > 3%) (|Smith et al.ll2002f l. Supporting this 
conclusion, Chandra and XMM-Newton observations of 
red 2MASS AGN find spectra that are eit her unusually 
hard or absorbed by Ng = 10^^ -10^^ cm~^ (jWilk es et al.l 
[200l[2005tlUrrutia et al"][2005l ). 

Unlike UV-excess based surveys, the Sloan Digital Sky 
Survey (SDSS) allows detection of red quasars in the op- 
tical band because the algorithm for selecting candidates 
for spectroscopy uses a four-dimension al multicolor se- 
lection criterion ( Richards et al.l l2003l) similar to tha t 
of the 2dF QSO Redshift Survev (jCroom et al.l[200l . 
This allows the SDSS to select all objects lying out- 
side of the stell a r locu s, including atypically red quasars. 
[Richards et al.l (j2003l ) explored the color distribution of 
SDSS quasars and found that while a population of in- 
trinsically red quasars exists, the majority of the red tail 
of the color distr ibution is explained by dust reddening. 
In a later study ([Hopkins et al.l[20M ). the curvatures of 
the SDSS spectra were found to be best-fit by a dust 
extinction cur ve similar to t hat of the Small Magellanic 
Cloud (SMC) (jPrevot et al .iri984). However, only mildly 
dust-reddened quasars can be detected (E(B-V) < 0.5) in 
the SDSS before dust extinction causes t hem to drop be- 
low t he detection threshold of the survey (jRichards et al.l 
[20nl . 

While dust-reddening is a common explanation for red 
colors in quasars, other possible explanations include: 
(1) host galaxy contamination, (2) absorption by the 
Lya forest, (3) optical synchrotron emission superim- 
posed on a normal, blue spectrum t o create a red spec- 
trum (Serjeant fc Rawlings l 119961: iFrancis et all 120001 : 
IWhiting et al.ll200ll l. and (4) an intrinsically red con- 
tinuum. The redshift selection in this paper minimizes 
contamination via options (1) and (2). 

X-ray observations used in conjunction with optical 
information can constrain the amount of intrinsic ab- 
sorption in a source, thereby allowing an investigation of 
dust-reddened versus intrinsically red optical continua. 
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Previous studies have found intrinsically red quasar can- 
did ates in addition to dust-reddened quasars. For exam- 
ple, |HaiOtIal] ()2006 l) selected a sample of 12 red quasars 
from the SDSS, but Chandra X-ray observations show ev- 
idence for NH absorption (Nh ~ 10^^ cm^^) for only 4 
of their 12 quasars. The other seven quasars show no ev- 
idence for even moderate intrinsic absorption and three 
of these show evidence for an intrinsically red optical 
continuum. 

In a similar vein, iRisaliti et all (|2003f ) observed 16 
optically-selected. X-ray weak quasars with Chandra. 
While the sample was not selected to be red, the aver- 
age color of the sample is redder than that of the parent 
Hamburg Quasar Survey {A{B-R) ^ 1 vs. A{B - R) - 
0.5). The X-ray weak quasars have a flat average photon 
index, < T >= 1.5 (where T = -a -I- 1 for oc ly") 
and 12 of the 16 quasars have steeper optical-to-X-ray 
indices than normal. (The optical-to-X-ray continuum is 
characterized by Uox = log(F2fcey/F™) / logiv2keV / J^uv) , 
where F2keV and F„„ are the intrir isic flux densities 
at 2 keV and 2500 A, respectively ()Tananbaum et aTl 
19791). 1 Whi le X -ray absorption remains a possibility, 
Risaliti et al.l ()2003f ) argue that these quasars are intrin- 
sically underluminous in the X-rays. The evidence, how- 
ever, is not conclusive due to the low signal-to-noise for 
many of the Chandra spectra. 

Larger surveys including high quality optical and X-ray 
spectra are now possible with the advent of the SDSS and 
the XMM-Newton archive of X-ray observations, both of 
which cover large areas of the sky. SDSS data can con- 
strain dust-reddening in the optical, while XMM-Newton 
observations constrain absorption in the X-rays. In this 
paper we investigate the optical and X-ray properties of 
17 red SDSS quasars. §2 introduces the data sources and 
sample selection and §3 covers the optical and X-ray re- 
duction and data analysis. In §4, we classify the quasars 
as dust-reddened, intrinsically red or unclassified, and in 
§5 we discuss low accretion rates as a possible explana- 
tion for the steep optical/UV spectra of intrinsically red 
quasars. We assume throughout the paper that Ho = 70 
km s-i Mpc-\ = 0.3, and f^A = 0.7. 

2. SAMPLE SELECTION 

2.1. Data Sources: SDSS & XMM Newton 

The third version of the SDSS quasar catalog 
(jSchneider et al.|[2005( ) is based on Data Release 3 (DR3), 
which covers an area of 5,282 deg^ (fAbazaiian et ahl 
[2005h . The quasar catalog contains 46,240 objects from 
the SDSS that have: (1) 15.0 < i < 22.2; (2) M, < -22.0; 
(3) at least one emission line with FWHM greater than 
1000 km s~^, or (4) unambiguous broad absorption line 
characteristics. SDSS photometry in the u, g, r, i and z 
bands covers 3,250 - 10,000 A. SDSS spectroscopy covers 
3,800 - 9,200 A (so the u band and part of the z band 
lie outside the spectral range), with a spectral resolving 
power of ~ 2,000. 

The SDSS DR3 quasar catalog overlaps with 1% of 
the area of archival XMM-Newton observations. XMM- 
Newton imaging is carried out by the three European 
Ph oton Imagin g Camer as (EPIC): MOS-1, MOS-2, and 
pn (jTurner et a l. 2001; Strueder et aI1l2001f ). The XMM- 
Newton telescope is good for retrieving serendipitous X- 
ray spectra of SDSS quasars due to a large collecting area 
(922 cm2 for MOS and 1,227 cm^ for PN at 1 keV), large 



field of view (33' x 33' for MOS and 27.5' x 27.5' for 
PN), and good spectral resolution (E/AE ~ 20 - 50 for 
both MOS and PN) The XMM-Newton spectral range 
covers the 0.5 - 12 keV band. 

The SDSS and XMM-Newton archives are well- 
matched in sensitivity. A source with a 2-10 keV flux 
of 10^^"* erg s^^ cm~^ and a standard radio quiet spec- 
trum (r ~ 2.0) produces ~ 100 counts in ~ 10 ksec, 
sufficient for a crude spectral fi t. For a typical value of 
= -1.5 (|Vignah et al.l |2003(), a 10"" ergs cm'^ s'^ 
quasar would have r = 19.37 while an X-ray quiet quasar 
with aox = —2.0 would have r = 15.63. Both of these 
magnitudes lie comfortably in the SDSS range. 

2.2. The SDSS-XMM Red Quasar Sample 

We define our red sample of [Schneider et al.l ()2005f ) 
quasars lying in XMM-Newton fields by selecting the 
most extreme red colors, g — r > 0.5. By further re- 
stricting our search to moderate redshifts (0.9 < z < 
2.1), we minimize host galaxy contribution at low red- 
shifts (z < 0.5) and Lya forest absorption at high red- 
shifts (z > 2.5). This selection results in a sample of 17 
q uasars. 

[Richards et al.l (|2003D . hereafter R03, defines dust- 
reddened quasars via their relative colors. Relative colors 
compare a quasar's measured colors with the median col- 
ors in its redshift bin, where redshift bin sizes are 0.1 in 
redshift, so that A{g — i) = [g — i) — < {g — i) >z. 
The use of relative colors corrects for the effect of typical 
emission lines on the photometry in a particular band. 
The relative colors of the SDSS quasars match a Gaus- 
sian distribution on the blue side but require the addition 
of a tail on the red side (Fig. 3 in R03). R03 identi- 
fies quasars belonging to the red tail as A(g — i) > 0.2, 
since the A{g — i) color correlates best with the pho- 
tometric spectral index. The A(g — i) color has been 
corrected for an SMC-like extinction curve with E(B - 
V) = 0.04, a value typical of high Galactic latitudes, for 
which F,{g - i) — 0.07. In practice, this yields a minimum 
of A{g — i) ~ 0.35 for the quasars with redshift z~1.5. 
The {g - i) and {g - r) colors behave similarly with red- 
shift (Fig. 2 of R03), and the maximum median color is 
0.2 for 1 < z < 2. In Figure [U we show a histogram of 
the relative (g - i) colors for our sample. Only two have 
A{g — i) < 0.35. Thus, the g - r > 0.5 color selection 
mimics A(g — i) > 0.35, though we do not use relative 
colors. 

Our red quasar sample is unbiased in X-ray loudness 
since we do not require X-ray detections for selection. 
Nevertheless, the X-ray detection rates are high: 94% 
for 2a detections (16/17) and 76% for > 3a detections 
(13/17), although most of the quasars have a relatively 
low X-ray signal-to-noise: 10/17 quasars have S/N < 10. 

Table [T] lists the sample, giving the source ID number 
used within this paper, the full SDSS name, radio clas- 
sification. X-ray S/N (which is averaged if multiple ob- 
servations are available), redshift. Galactic Nh, and the 
absolute i-band magnitude M^. A source is radio lo ud if 
F5GHZ I F4400 = Rl > 10 (|Kellerman et al.lll989D . To 
calculate -Rl, a power-law is interpolated between the 
optical magnitudes to get Fa(4400), and the 1.4 GHz 
radio flux is obtained from FIRST survey detections 

* http: / /imagine. gsfc.nasa.gov/docs/sats_n_data/missions/xmm. html 
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Fig. 1. — The A{g — i) color distribution of the red sample. The 
majority of sources in this sample meet the [Richards et al.l (12003 ) 
requirement that A(g — i) > 0.35 for dust-reddened quasars. 



(jWhite et al.1119971 ). which are extrapolated to 5 GHz 
using a radio power-law an = -0.8. All the quasars lie 
in the area covered by the FIRST survey, so if there is 
no detection, we use the 5a upper-limit on the 1.4GHz 
radio flux to extrapolate to 5 GHz. 

For five of 17 sources, there are two X-ray observations. 
In these cases, both observations were processed and used 
in our analysis. 

2.3. SDSS Data 

SDSS photometry was extracted from the online 
database^ via SQL-based queries in CasJobs. The in- 
dividual, calibrated SDSS spectra were also downloaded 
from the SDSS database. Table [5] gives the source ID, 
the abbreviated SDSS name, the g — r color (which 
determined selection), the apparent i band magnitude, 
and the full- width half-maximum (FWHM [km s^^]) of 
the broad Mgll component. The Mgll information is 
taken from lShen et al. I ()2007[ ) for all but three objects. 

(12007D simultaneously fit the Fell emission, 
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a continuum power-law, and broad and narrow Gaus- 
sians to the Mgll line. The spectral signal-to-noise is too 
low for two objects (#13 and #15) to use this method. 
Another ob ject (#17) has goo d S/N, but line values are 
not given in lshen et aFI ()2007[) . We fit the Mgll lines for 
these three sources by estimating the continuum values 
by eye and then simultaneously fitting broad and narrow 
Gaussians with IRAF^. 

There is one broad-absorption line quasar (BAL) in 
the sample, SDSS J0944-h0410 (#1), and two narrow 
absorption-line (NAL) quasars, SDSS J1226-t-0130 (#8) 
and SDSS J1435-)-4851 (#9). Only five of the 17 quasars 
have the CIV line, which is the most reliable marker of 
BAL/NAL presence, visible in their spectra. Mgll, the 
BEL seen in all 17 quasars, rarely shows either BALs or 
NALs. So the true number of BALs and NALs in the 
sample could be even greater. 

2.4. X-ray Data Reduction 

Using the XMM-Newton Science Analysis System, 
SAS v7.02^, we filtered the XMM observations to remove 
time intervals of flaring high-energy background events 

^ http:/ /www. sdss.org/dr3/access/index. html 
® http://iraf.noao.edu/ 
http://xmm.csac.csa.int/sas 



using the standard cut-off of 0.35 cts s~^ for the MOS 
cameras and 1.0 ct s^^ for the FN camera. The event 
rate plots were visually examined to determine that these 
cut-offs were appropriate. We then extracted source and 
background regions for spectral analysis. The SAS task 
eregionanalyse was used to select the source radius with 
the optimal signal-to-noise, with typical radii ranging 
from 10" to 60" . Background regions were defined by eye, 
avoiding obvious X-ray sources and chip edges. These 
regions were typically a circle of radius 2000-2500 pixels 
(100-125"), selected to lie on the same chip as the source 
and as close to the source as possible without overlapping 
the source extraction region. 

Where possible, observations were retrieved for all 
three XMM EPIC CCDs. In seven observations, the 
source lies in a bad region in one of the cameras, either in 
a strip between two chips or, because the MOS and FN 
cameras have different shapes, the source may lie outside 
the field of view in one of the cameras. In an eighth case, 
the source is within the field of view of all three cameras, 
but in the MOS-2 camera, most of the events are filtered 
out while removing the flaring high-energy background. 
In all these cases (7 of 22 observations), we use the re- 
maining images from the other cameras for analysis. 

Table [3] contains data for the XMM observations, or- 
dered by increasing S/N. The table includes the source 
ID, XMM identification number, observation date, net 
exposure time, net counts. X-ray S/N and the cameras 
available for analysis. 

3. ANALYSIS 

3.1. Optical Spectra and Continuum Dust Reddening 

Fits 

We wrote an IDL procedure to fit dust-reddening to the 
optical continuum. The procedure calls a de-reddening 
routine, FM.UNRED.FRO (Landsman 1999), which uses 
the Fitzpatrick (1999) parameterization to characterize 
the optical and UV extinction curve. 

Each source spectrum is first shifted to the rest-frame. 
A composite quasar spectr um created from SDSS quasars 
(jVanden Berk et "all 1200 1^ is then redde ned incremen- 
tally with the SMC extinction curve ijGordon et al.l 
l2003f l and normalized to the red quasar spectrum by 
averaging the cont i nuum between 3020 and 3066 A 
( Van den Berk et al] l2001h in order to calculate the 
value. The best-fit parameters are those for which 
reached a minimum. The total was allowed to vary 
by 2.7 to get the 90% confidence interval errors. 

The SMC, rather than Galactic, extinction curve is 
used because of the absence of the A2175 feature in 
AGN spectra ("Fitman et aTl 120 00: Kuh n et al."l l2001l : 
ICzernv et all [2004: Ho pkins et al. 2004). The optical 
dust-reddening (hereafter denoted as E(B- V)spec) is al- 
lowed to vary from E(B- V)spec = 0.0 to E(B- V)spec = 
1.0. Simulations have shown that quasars with E(B-V) 
> 0.5 are not expecte d from the SDSS selection process 
([Richards et al.l [2003fl. so the fits to the spectral contin- 
uum cover the full E(B-V) parameter space. 

The Small Blue Bump, a blend of Balmer continuum 
emission and Fell lines, extends from 2000 - 4000 A 
([Wills. Netzer fc Wills1ll985l ). Our color selection crite- 
ria {g — r > 0.5) bias our sample towards sources with 
atypically large Fell emission between ~ 2200 - 2700 A. 
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TABLE 1 

SDSS/XMM-Newton Red Quasar Sample 



Source 
ID 


SDSS name 


Rl" 


(S/N),'' 


z 




Mi 


1 


SDSS J094440.43+041055.5 


RQ 


0.17 


1.98 


3.63 


-25.3 


2 


SDSS J023217.68-073119.9 


RQ 


3.23 


1.16 


3.05 


-23.6 


3 


SDSS J165245. 77+394722.3 


RQ 


1.89 


1.19 


1.64 


-24.2 


4 


SDSS J072843. 03+370835.0 


RL (Ri = 34.8) 


2.87 


1.40 


6.00 


-24.0 


5 


SDSS J100201. 50+020329.4 


RQ 


2.85 


2.01 


2.59 


-25.0 


6 


SDSS J221723. 29-082139.4 


RQ 


2.93 


1.25 


5.30 


-24.2 


7 


SDSS J153322. 80+324351.4 


RQ 


4.04 


1.90 


2.05 


-24.8 


8 


SDSS J122637.01+013016.1 


RQ 


4.69 


1.55 


1.84 


-25.2 


9 


SDSS J143513. 89+484149.3 


RQ 


4.46 


1.89 


2.07 


-24.7 


10 


SDSS J221719. 55-081234.4 


RQ 


5.88 


1.46 


5.31 


-24.1 


11 


SDSS J111452. 84+531531. 7 


RQ 


6.33 


1.21 


0.923 


-24.1 


12 


SDSS J113342. 73+490025.8 


RQ 


8.65 


1.30 


1.59 


-24.6 


13 


SDSS J125456. 82+564941.4 


(Ri <18.3) 


14.43 


1.27 


1.26 


-22.8 


14 


SDSS J095918. 70+020951.4 


RL (Rl = 36.3) 


21.87 


1.16 


2.64 


-23.4 


15 


SDSS J095857.34+021314.4 


RQ 


27.94 


1.02 


2.64 


-23.0 


16 


SDSS J105316. 76+573550.7 


RQ 


35.21 


1.21 


0.556 


-23.9 


17 


SDSS J091301. 01+525928.9 


RQ {Rl = 4.9) 


43.48 


1.38 


1.60 


-26.9 



TABLE 2 
Red Quasar Sample in the SDSS 



ID 


SDSS name 


g _ r 


1 


Mgll FWHM 
(km s-i) 


1 


SDSS J0944+0410 


0.50 


18.22 


8,800 


2 


SDSS J0232-0731 


0.53 


19.25 


11,700 


3 


SDSS J1652+3947 


0.54 


18.65 


8,800 


4 


SDSS J0728+3708 


0.74 


19.09 


9,700 


5 


SDSS J1002+0203 


0.55 


18.58 


10,800 


6 


SDSS J2217-0821 


0.74 


18.69 


7,600 


7 


SDSS J1533+3243 


0.97 


18.69 


6,700 


8 


SDSS J1226+0130 


0.59 


18.01 


7,700 


9 


SDSS J1435+4841 


0.50 


18.78 


7,000 


10 


SDSS J2217-0812 


0.62 


19.06 


5,800 


11 


SDSS J1114+5315 


0.76 


18.77 


6,400 


12 


SDSS J1133+4900 


0.52 


18.42 


9,100 


13 


SDSS J1254+5649'' 


0.57 


20.12 


8,200 


14 


SDSS J0959+0209 


0.60 


19.40 


23,200 


15 


SDSS J0958+0213'' 


0.59 


19.66 


10,500 


16 


SDSS J1053+5735 


0.55 


18.95 


5,400 


17 


SDSS ^913+5259=^ 


0.64 


16.20 


6,000 



Note. — FWHM measureme nts for the Mgll broad emission 
lines are from lShen et al. I II2007I) . 

Five of 17 sources have atypical Fell emission. Excluding 
the region from 2200 - 3000 A from the chi-square calcu- 
lation for all the sources accounts for both Fell emission 
and the 2798 A Mgll line, so that the reddening fit to the 
continuum is not biased by unusually strong line emis- 
sion. The reddening fit was not sensitive to other broad 
emission lines, such as the 1549 A CIV line, so we did not 
omit any additional lines from the fit. The BAL quasar 
(SDSS J0944+0410, #1) was not fit with reddening be- 
cause numerous absorption lines and a strong Fe II bump 
left no reliable continuum points with which to calculate 

We also fit a power-law (Fi, oc to the optical 

continuum. All major emission lines are removed before 
performing the fit, as is the region dominated by the 
blended Fell and Balmer emission. A power-law will not 
be a good fit to spectral continua with obvious dust cur- 
vature, but for (E(B-V)spec < 0.1), dust-reddening will 
look similar to a red power-law. Section §3.2 discusses 
this issue in more detail. 

Figure [2] plots the SDSS source spectra in Fa (10^^'' 

ergs cm" -2 s-1 A ) VS. A (A) space. Overplotted are 



the original FVanden Berk et all (|2001h template (dotted 
blue line), the reddened template (dashed red line), and 
the optical power-law (solid green line) . Table [4] gives 
the reddening results: E(B-V)spec, the expected intrinsic 
gas column density {^H.o vt) obtained by a.pplying the 
Gala ctic gas-to-dust ratio (jBohlin et al.lll978l : lKent et al.l 
|1991| ) , the intrinsic gas column density from the X-ray fits 
(N_f/,£c, see §3.3), and the optical power-law index, aopt- 

Reddening fits to the optical spectra require E(B- 
V)spec in 14 out of 17 sources at> 3tT. For these 14 
sources, E(B-V)spec ranges from 0.04 to 0.25, with a 
weighted mean of < E{B — V)spec > = 0.1 correspond- 
ing to < NH,opt > = 5.8 X 10^° cm~2, using the Galactic 
gas-to-dust ratio. 

3.2. Testing the Optical Continuum for Weak Curvature 
with Relative Colors 

We find that reddening is a good fit to most of the 
optical spectra. However, for the redshifts of this sample 
(z > 1), the optical continuum will not have strong dust- 
induced curvature for mild reddening (E(B-V) spec < 0.1). 
Under these conditions, a dust-reddened optical contin- 
uum will look similar to a red power-law. Figure[2]in this 
paper shows the similarity between mild dust-reddening 
(dashed red line) and a red power-law (solid green line). 
The curvature in the relative broad-band photometry can 
help differentiate these two possibilities because the u 
band (for which the FWHM ranges from 3251 to 3851 A) 
is mostly outside the spectral range (3800 - 9200 A) . The 
u band is a factor 3.6 more absorbed by dust (for E(B-V) 
= 0.1) than the r band, where the spectra are centered. 

Relati ve colors give an idea of the shape of the con- 
tinuum (|Hall et al.ll2006l) . A quasar with a typical g — i 
color, for example, will have lS.{g — i) = 0, while a quasar 
with red [g — i) color will have A{g ~ i) > 0. A dust- 
reddened quasar with typical emission lines will have 
A{u — r) > A{g — i) > A(r — z), because curvature 
increases on the blue side of the continuum. A typical 
blue power-law continuum, on the other hand, will have 
A{u — r) ~ A{g — i) ~ A{r — z). In the case of a power-law 
continuum that is redder than average, the relative colors 
will go as A{u ~r) < A(g — i) > A(r — z). Therefore, an 
object with A{u—r) — A{g~i) < displays dust-reddened 
curvature, while an object with A{u — r) — A{g — i) > 



The Reddest DR3 SDSS/XMM-Newton Quasars 



5 



TABLE 3 
XMM-Newton Observation Log 



Source 


XMM ID 


Date 


Net Exp Time'' 


Net Counts'' 


S/N 


Cameras 


ID 






(ks) 


(counts) 


PN Ml M2 


1 


0201290301 


2005-6-09 


56.7 


0.03 


0.2 


-XX 


2 


0200730401 


2005-1-30 


122.8 


40.3 


3.2 


XXX 


3a 


0113060401 


2002-07-14 


4.55 


7.3 


1.2 


-XX 


3b 


0113060201 


2003-9-03 


5.0 


19.2 


2.6 


-XX 


4 


0145450501 


2003-12-04 


9.2 


25.1 


2.9 


XXX 


5a 


0203360401 


2005-12-10 


29.6 


39.0 


2.8 


XXX 


5b 


0203360801 


2005-12-07 


89.4 


11.8 


2.9 


XX- 


6 


0009650201 


2002-9-20 


40.2 


31.2 


2.9 


-XX 


7 


0039140101 


2003-8-29 


25.3 


70.8 


4.0 


XXX 


8 


0110990201 


2002-8-11 


32.3 


70.3 


4.7 


XXX 


9a 


0110930401 


2004-2-04 


21.1 


77.5 


4.2 


XXX 


9b 


0110930901 


2004-2-04 


17.7 


98.3 


4.7 


XXX 


10 


0009650201 


2002-9-20 


40.2 


84.5 


5.9 


-XX 


11 


0143650901 


2004-5-20 


17.1 


135.3 


6.3 


XXX 


12 


0149900201 


2004-12-15 


36.0 


196.3 


8.7 


-XX 


13 


0081340201 


2002-11-02 


60.1 


750.0 


14.4 


XXX 


14 


0203361801 


2005-1-16 


79.3 


1486.8 


21.9 


XXX 


15 


0203361801 


2005-1-16 


79. 


2696.2 


27.9 


XXX 


16a 


0147511701 


2004-2-01 


282.4 


4185.3 


33.7 


XXX 


16b 


0147511801 


2004-2-01 


282.2 


4775.3 


36.8 


XXX 


17a 


0143150301 


2004-9-06 


20.8 


3687.3 


33.6 


XXX 


17b 


0143150601 


2004-9-06 


47.9 


9130.1 


53.4 


XXX 



TABLE 4 
Optical Reddening Fits 



ID 


E(B-V)spec 


(1020 cm-2) 


(1020 cm-2) 




1 








-1.49 


2 


0.03±0.02 


1.7±1.2 




-0.92 


3 


0.08±0.02 


4.6±1.2 




-1.05 


4 


0.25±0.02 


14.5±1.2 




-2.92 


5 


0.052±0.005 


3.0±0.3 




-1.36 


6 


0.17±0.02 


9.9±1.2 




-1.95 


7 


0.12±0.02 


7.0±1.2 


<696 


-1.85 


8 


0.15±0.01 


8.7±0.6 


<1340 


-1.79 


9 


0.13±0.01 


7.5±0.6 


<677 


-2.07 


10 


0.11±0.02 


6.4±1.2 


<184 


-1.68 


11 


0.18±0.03 


10.4±1.7 


<89.0 


-1.46 


12 


0.04±0.01 


2.3±0.6 


<76.3 


-0.87 


13 




t: + 6.4 
°-*-4.6 


<23.1 


-1.44 


14 


0.12±0.03 


7.0±1.7 


<5.49 


-1.05 


15 


16+009 


n 0-I-5.2 
M.O_4 1 


<5.39 


-1.94 


16 


0.11±0.02 


6.4±1.2 


<4.7 


-1.17 


17 


n n8fi+0-004 




<2.9 


-1.31 



displays a red power-law. We use this criterion to clas- 
sify quasars as intrinsically red, as shown in Figure [S] 
(discussed in §4.2). 

Since relative colors compare the observed quasar col- 
ors to the mean colors for quasars at that redshift, atypi- 
cal emission lines will affect the relative color results. For 
the redshift range 1 < z < 2, the most significant lines 
are Mgll and the Fell emission line blend. To correct for 
atypical Mgll lines, we calculate the equivalent width of 
Mgll from the spectrum. From this, we subtract the av- 
erage Mgll equivalent width obtained from the Vanden 
Berk composite spectrum (2001). The residual Mgll line 
flux is calculated and added or subtracted from the band 
where Mgll is found. To correct for the Fell emission 
line blend, we redden the composite spectrum, shift it 
to the observed frame, and normalize to the source spec- 
trum. We then subtract this spectrum from the source 
for the range 2200-2700 A and calculate the equivalent 
width of the residual Fell line. We add or subtract the 
residual flux from the bands where Fell is found. Since 



the u band is outside of the spectral range, we cannot 
correct for atypical emission lines in this band. The rel- 
ative colors and the optical continuum shapes are listed 
in Tabled 

3.3. X-ray Spectral Fits 

We made fits to the extracted spectra using the Sherpa 
package^ within CIAO^. For each source, the available 
MOS-f PN spectra were fit simultaneously by linking pa- 
rameters. We used a conservative spectral range of 0.5 - 
10 keV for fitting. The observations were fit according to 
their S/N, with more complicated models (models A-E) 
being applied as S/N increased. Table [6] summarizes the 
model characteristics. All models applied in this paper 
include local absorption fixed to the Galactic hydrogen 
column density {NH,gai) for a source's coordinates. Val- 
ues for NH,gai were taken from WebPIMMS^". 

For six observations with S/N < 4, we fit Model A, 
where we freeze the photon index of the X-ray power- 
law to F = 1.9. We leave normalization free to vary in 
order to get the flux or a 90% upper-limit on the flux. 
For five observations with 2 < S/N < 4, we fit Model B, 
a single power-law with free F and normalization. Mod- 
els A and B use the Cash statistic for fitting, which we 
discuss further below. 

For 14 observations with S/N > 4, we fit Model C, a 
single power-law, but using Chi-square statistics rather 
than Cash. For the same observations, we also fit the 
spectra with Model D, a power-law with intrinsic ab- 
sorption. Finally, for five observations with S/N > 25, 
we fit model E, a power-law with intrinsic absorption, 
plus a Gaussian for the Fe Ka line. The Gaussian line 
energy is set to 6.4 keV shifted to the quasar's reference 
frame. For sources with S/N > 4 (where models C, D 
and E are applied), the Chi-square statistic is used and 
the data are binned by at least 15 counts/ bin. 

* http://cxc.harvard.edu/sherpa/threads /index. html 
9 http://cxc.harvard.edu/ciao/ 

-"^0 http:/ /heasarc. gsfc.nasa.gov/Tools/w3pimms. html 
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TABLE 5 

Relative Colors and the Optical Continuum Shape 



ID 


A( 


u-r) 




Al 


IS-i) 


A(r-z) 


Continuum Shape 






1 


0.97 


±0.06 


0, 


,68 


±0.06 


0.13 ±0.04 


Dust 






2 


0.96 


±0.07 


0, 


,25 


±0.04 


0.03 ±0.07 


Dust 






3 


0.79 


±0.04 


0, 


,41 


±0.03 


0.29 ±0.04 


Dust 






4 


1.26 


±0.08 


0. 


,81 


±0.03 


0.65 ±0.04 


Dust 






5 


0.63 


±0.04 


0. 


,71 


±0.04 


0.16 ±0.04 


Red power-law (oopt = 


-1 


.36) 


6 


1.06 


±0.07 


0, 


,83 


±0.03 


0.52 ±0.05 


Dust 






7 


1.10 


±0.06 


0, 


,69 


±0.02 


0.23 ±0.04 


Dust 






8 


1.13 


±0.06 


0. 


,74 


±0.04 


0.58 ±0.04 


Dust 






9 


0.99 


±007 


0, 


,77 


±0.03 


0.56 ±0.04 


Dust 






10 


0.77 


±0.07 


0, 


,74 


±0.03 


0.41 ±0.06 


Undefined 






11 


0.80 


±0.07 


0. 


,74 


±0.05 


0.47 ±0.04 


Undefined 






12 


0.14 


±0.04 


0, 


,24 


±0.03 


0.06 ±0.04 


Red power-law («opt = 


-0 


.87) 


13 


-0.07 


±0.07 


0, 


,54 


±0.05 


0.56 ±0.1 


Red power-law («opt = 


-1 


.44) 


14 


0.60 


±0.07 


0, 


,54 


±0.03 


0.38 ±0.08 


Undefined 




15 


0.54 


±0.08 


0, 


,63 


±0.04 


0.40 ±0.09 


Red power-law (aopt = 


-1 


.94) 


16 


0.33 


±0.04 


0, 


,44 


±0.03 


0.26 ±0.04 


Red power-law {aopt = 


-1 


.17) 


17 


0.25 


±0.02 


0, 


,40 


±0.02 


0.20 ±0.03 


Red power-law {aopt = 


-1 


.31) 



Note. — The relation between the relative colors gives the optical continuum shape from the photometry. A(u - r) > A(g - i)> A(r - 
z) indicates dust, but A(u - r)< A(g - i)> A(r - z) indicates a red power-law. 



TABLE 6 

S/N-DEPENDENT MODELS APPLIED TO XMM OBSERVATIONS 





Model A 


Model B 


Model C 


Model D 


Model E 


S/N 


S/N < 4 


2 < S/N < 4 


9/N > 4 


9/N > 4 


9/N > 25 


Model 


Cash fit 


Cash fit 


fit 


fit 


fit 


Description 


Fixed power-law 


Free power-law 


Free power-law 


Free power-law 


Free power-law 




(r = 1.9) 












Gal. abs. 


Gal. abs. 


Gal. abs. 


Gal. abs. 


Gal. abs. 




(fixed) 


(fixed) 


(fixed) 


(fixed) 


(fixed) 










Intrinsic 


Intrinsic 










abs. (free) 


abs. (free) 










Fe Ko line 



Note. — A total of 22 observations cover 17 sources. 

For S/N < 4, Chi-square is not an appropriate statis- 
tic because there are not enough counts per bin, so both 
models A and B are fit using the more time-consuming 
Cash statistic (1979) which gives more reliable results 
for low-count sources (Nousek & Shuc 1989). When us- 
ing Cash statistics, the source and background counts 
are binned by one count/bin, to ensure that there are 
no empty bins. The background is not subtracted and 
is instead fit simultaneously with the source. The XMM 
background is fit with the three components described 
in the XMM Users Handbook^ ^: a power-law fixed to 
r = 1.4 plus Galactic absorption (to account for the ex- 
tragalactic X-ray spectrum), a broken power-law, with 
the break energy fixed to 3.2 keV (to account for the 
quiescent soft proton spectrum), and lines at 1.5 keV 
for the MOS cameras and at 1.5 and 8 keV for the PN 
camera (to account for cosmic-ray interactions with the 
detector). Since this background model has five free pa- 
rameters for each camera, constraining all fit parameters 
results in large errors in the best-fit parameters. 

Table [7] lists the source ID, the model used for each 
source, XMM observation ID, the best-fit spectral slopes, 
flux values from 0.5-2 keV and 2-10 keV, the rest-frame 2- 
10 keV luminosity, and the or Likelihood of Fit values 
and degrees of freedom for the best fit. For sources with 
S/N<4, Table[7]hsts the best- fit parameters for Model 
A, where the spectral slope is fixed to 1.9, leaving the 



normalization free to vary. (For sources with S/N < 2, 
the best-fit upper-limits are listed.) We list the best-fit 
spectral slopes from Model B in Table [5] for reference, 
but due to the large spectral slope errors we do not use 
these values for further analysis. 

For sources with S/N > 4, we determine the best- fit 
model via the F-test^^, which measures the significance 
of the change in chi-square as components are added to 
a model. In all cases, the more complex models did not 
have a significantly lower chi-square value, so for S/N > 
4, Table [7] lists the parameters for the power-law model 
(model C). 

Tables [5] and [Sllist the best-fit spectral slope and the 
upper-limits on other spectral parameters for rejected 
models B and D (Table [8|) and E (Table [9]). Figure [3] 
shows the data, fits, residuals and la, 2a and 3a F-N^f 
contours for the sources with S /N > 4 fitted with model 
D (power- law -I- intrinsic absorption). 

3.4. X-ray Absorption 

Of the 11 sources (and 14 observations) with S/N > 4, 
none show significant intrinsic absorption, and the F-test 
demonstrates that all sources prefer a simple power-law 
model. Nh upper limits (Model D) for the five highest 
S/N sources range from 3 x 10^° cm^^ to 2 x 10^^ cm~^ 
at the 90% confidence level. 



Absorption can also be indi cated when th e best- fitting 
power-law is unusually hard. iMateos et al.l (|2005D finds 

http: // xmm. vilspa. esa.es/ external / xmm_user .support / documentation 

http:/ /cxc. harvard.edu/ciao/ahclp/ftest. html 
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TABLE 7 
Best-fit X-ray spectral fits 



Source 
ID 


Model 


XMM ID 


r 


Po.5-2kev'^ 




L2—10keV*' 




1 


A 


0201290301 


1.9" 


<0.024 


<0.034 


<9.5 X 10*^ 


1791/1599 


2 


A 


0200730401 


1.9 


0.24 


0.36 


5.2 X 10^2 


8904/5766 


3a 


A 

A 


011dUDU4Ul 


1.9 


<0.51 


<3.72 


<.3.D X lU 


lz2o/loUo 








1 Q 


2 25 


^ 1 n 

O. J-U 


4 8 X lO^^ 


1 ni s /I fil 4 


4 


A 


0145450501 


1.9 


0.016 


13.4 


5.2 X 10"^ 


1608/1619 


5a 


A 


0203360401 


1.9 


0.32 


0.46 


1.3 X 10*3 


17119/5749 


5b 


A 


0203360801 


1.9 


0.53 


0.74 


2.1 X lO''^ 


10447/4946 


6 


A 


0009650201 


1.9 


0.52 


2.96 


3.2 X lO''^ 


1831/1636 


7 


C 


0039140101 


l.8tj];t 


1.95 


3.05 


7.7 X lO'*''* 


4.3/6 


8 


C 


0110990201 




1.36 


7.49 


7.0 X lO'*^ 


9.7/7 


9a 


C 


0110930401 


^•^-0.4 


2.78 


7.96 


1.3 X lO"** 


2.6/5 


9b 


c 


0110930901 




4.57 


27.2 


2.7 X 10** 


3.7/7 


10 


c 


0009650201 


1 O+0-4 


3.05 


4.90 


9.4 X lO*'"* 


2.6/6 


11 


c 


0143650901 


-, r+0.3 
-"-•"-(). 3 


5.03 


13.3 


1.5 X lO** 


5.9/9 


12 


c 


0149900201 




8.58 


16.2 


3.0 X 10*4 


8.1/15 


13 


c 


0081340201 


1 6+°-i 


9.59 


20.7 


2.8 X 10** 


27/40 


14 


c 


0203361801 


-1 7c:+0.09 


6.52 


11.9 


1.6 X 10** 


54/63 


15 


c 


0203361801 




17.0 


29.5 


3.4 X 10** 


75/107 


16a 


c 


0147511701 


-, 77+0.06 
' ' -0.05 


12.1 


20.1 


2.9 X 10** 


100/165 


16b 


c 


0147511801 


1 79+0.06 
^- '^-0.05 


13.2 


23.5 


3.2 X 10** 


155/189 


17a 


c 


0143150301 


1 68+"-°5 


60.6 


120 


1.8 X 10*5 


94/142 


17b 


c 


0143150601 


1 79+0.04 
-^•'^-0.03 


57.9 


107 


1.7 X 10*5 


203/311 



TABLE 8 

Alternative X-ray Spectral Fits: r (S/N < 4) and N^^^. (S/N > 4) 



Model 


Source 


XMM ID 


r 








ID 






(1022 cm-2) 


B 


2 


0200730401 






8949/5767 


B 


3a 


0113060401 






1015/1614 


B 


4 


0145450501 


1 9+2.6 




1418/1619 


B 


5a 


0203360401 






17008/5750 


B 


5b 


0203360801 


9 c + 18.5 




11003/4946 


B 


6 


0009650201 


1 2+2'^ 




1822/1636 


D 


7 


0039140101 


-'-•''-0.4 


<7.0 


3.5/5 


D 


8 


0110990201 


^•"-0.4 


<13.4 


8.0/5 


D 


9a 


0110930401 


1 4+°-* 


<6.8 


2.6/4 


D 


9b 


0110930901 


1 0+"-* 


<10.2 


3.6/6 


D 


10 


0009650201 


1 9+"'* 


<3.3 


2.2/5 


D 


11 


0143650901 


-^-"-0.3 


<0.9 


5.9/8 


D 


12 


0149900201 


1 4+0-5 


<0.8 


7.3/14 


D 


13 


0081340201 


1 6+°-l 
^•"-0.1 


<0.2 


27/39 


D 


14 


0203361801 


1 7r-+0.09 
^- "^-0.09 


<0.05 


54/62 


D 


15 


0203361801 


1 7O+0.07 
^■'=-0.06 


<0.05 


75/106 


D 


16a 


0147511701 


1 77+0.06 
^- ' ' -0.05 


<0.06 


100/164 


D 


16b 


0147511801 


-I 79+0.06 
^•'^-0.05 


<0.04 


155/188 


D 


17a 


0143150301 


1 fiS+O-OS 


<0.04 


94/141 


D 


17b 


0143150601 


-1 79+0.04 

-^•'"'-o.os 


<0.02 


203/310 



Note. — For sources with S/N<4, model parameters from the Model B (Cash power-law) fits are given. For sources with S/N>4, we 
list 90% confidence upper-limits on Model D, power-law -|- intrinsic absorption. 
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TABLE 9 

Alternative X-ray Spectral Fits - Fe Ka line (S /N > 25) 



Source 


Model 


XMM ID 


Fe Ka FWRM* 


Fe Ka line flux^ 




ID 






(eV) 


(photons cm~^ 


15 


E 


0203361801 


<10 


<2.6e-6 


72.3/102 


16a 


E 


0147511701 


<220 


<3.6c-7 


100.2/160 


16b 


E 


0147511801 


<440 


<2.5c-6 


145.8/184 


17a 


E 


0143150301 


<190 


<3.8e-4 


92.2/137 


17b 


E 


0143150601 


<0 


<3.3e-6 


194.5/306 
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2. SD£5 J0252-073I E^., = 0.03, 





3. SPSS J1E5243947 Eg.y = O.QB, o„,= -1.05, z =1.19 



SPSS J072S + 370a = Q 25 Q„p,= -2.92, z =1.40 





6. SPSS J22l7-0g2l Eb_. = 



= -1.95, z =1,24 



0. SPSS J221 7-0812 E^.., = 0.1 1, o„,= -1,68, z =1.46 





Unclassified _ 








Fig. 2. — The SDSS source spectra are plotted as solid black 
line s, F\ (10"^''^ e rg s cm ~^ s~^ A) vs. A (A). The original 
IVandcn Berk ct al. ( 2001!) template (dotted, blue line) and the red- 
dened template (dashed, red line) are overplotted. Note that the 
BAL quasar (SDSS J0944-I-0410, #1) had no reliable continuum 
to fit reddening and so only the original template is plotted. The 
solid line beneath the spectrum marks the Fell+Mgll region ig- 
nored by the continuum-fitting program as described in the text. 
A power-law fit to the spectrum, omitting major emission lines, is 
also shown (solid, green line). 



an average X-ray photon index < F > 1.96 and in- 
trinsic dispersion crp = 0.4 for a sample of 1137 AGN 
found in XMM-Newton fields, with fluxes ranging from 
10-15 - 10-12 ergs cm-2 s-i. iMateos et all (120051 ) find 
a systematic hardening of the average X-ray spectra to- 
wards fainter fluxes, which they interpret as a higher 
degree of photoelectric absorption among fainter AGNs. 



TABLE 10 

X-RAY Strength Relative to the Optical 



ID 


J- uv 


F t> 

J- uv,corr 




■p d 

^ 2ke'V,corr 




f^ox ,corr^ 


1 


3.87 


6.04 


<2.45 


<2.45 


<-1.61 


<-1.69 


2 


3.19 


3.99 


0.33 


0.33 


-1.91 


-1.95 


3a 


4.39 


7.54 


<0.17 


<0.17 


<-2.07 


<-2.16 


3b 


4.39 


7.54 


3.08 


3.08 


-1.59 


-1.68 


4 


2.85 


15.05 


2.81 


2.81 


-1.54 


-1.82 


5a 


2.69 


3.83 


0.38 


0.38 


-1.86 


-1.92 


5b 


2.69 


3.83 


0.67 


0.67 


-1.77 


-1.83 


6 


5.00 


15.66 


0.15 


0.15 


-2.12 


-2.31 


7 


2.53 


5.82 


2.75 


29.0 


-1.52 


-1.27 


8 


5.65 


15.34 


0.77 


7.24 


-1.87 


-1.66 


9a 


2.10 


5.06 


1.87 


45.1 


-1.55 


-1.17 


9b 


2.10 


5.06 


1.38 


3.67 


-1.61 


-1.59 


10 


2.64 


5.67 


4.32 


10.0 


-1.45 


-1.44 


11 


5.34 


17.66 


4.57 


8.44 


-1.56 


-1.66 


12 


7.39 


9.88 


7.36 


10.6 


-1.54 


-1.52 


13 


1.33 


2.64 


9.72 


10.9 


-1.20 


-1.30 


14 


2.20 


5.09 


7.79 


8.02 


-1.32 


-1.46 


15 


2.21 


6.47 


21.9 


22.6 


-1.15 


-1.33 


16a 


4.42 


9.48 


14.7 


15.1 


-1.33 


-1.45 


16b 


4.42 


9.48 


15.2 


15.6 


-1.33 


-1.45 


17a 


53.67 


95.96 


63.5 


64.5 


-1.51 


-1.60 


17b 


53.67 


95.96 


63.8 


64.4 


-1.51 


-1.60 



Therefore, we define an unusually hard spectrum to have 
r < 1.5 at 90% confidence. Two of 14 sources with S/N 
> 4 have unusually hard spectral slopes by this mea- 
sure: SDSS J1226+0130 (#8) and SDSS J1435-H4841 
(#9). AU of the quasars with high enoug h S/N to fit T 
have b est-fit spectral slopes flatter than the IMateos et al.l 
(|2005l ) average. 

For the five sources with 2 < S/N < 4 (fit with model 
A), the S/N is not high enough to fit a power-law or 
intrinsic absorption. One indicator of X-ray absorption 
for these sources is X-ray weakness relative to the opti- 
cal luminosity, defined as aox < —1.8, where Uox is the 
spectral index from 2500 A to 2 keV. Such small values of 
aox are a common characteristic of BAL quasars, which 
are known to be heavily absorbed ob.jects in the X-rays 
(iMathur et al.l[T995HC^en fc Mathur Ill99()f) . 

To calculate aox for the red sample, we use the 2 keV 
flux (F2kev) given by the X-ray fits, and we derive the 
2500 A flux (Fuy) by interpolating a power-law between 
the nearest two of the ugriz magnitudes. We then calcu- 
late aox,corr by correcting aox for ab sorption. We dered- 
den F„„ using E{B - V)spec and the lPrevot et al.l ()1984D 
SMC extinction curve, where Ax = E(B-V) 1.39 A-^-^. 
For objects with X-ray S/N > 4, we do not find signifi- 
cant intrinsic absorption, so we use the 90% upper-limit 
on Nh to correct F2keV- For objects with X-ray S/N 
< 4, we cannot fit intrinsic absorption, so we make no 
correction to F2keV ■ Table [10] gives the measured and 
corrected F„„ and F2keV, olox (uncorrected), and aox,corT 
(corrected for absorption). 

Figure |4] plots aox against the log of the 2500 A lu- 
minosity (/„„) for the red sample. The unfilled circles 
give the measured aox and luv values for the red sample. 
The filled circles give aox.corr and luv.corr, the 2500 A 
luminosity corrected for absorption. Solid lines connect 
the uncorrected values to the corrected values for each 
source. The five sources with X-ra y S/N < 4 are marked 
in red. The iStrateva et al.l ()2005f ) sample and aox — luv 
correlation are plotted as black points and a solid line, 
respectively, for comparison. The dotted line marks aox 
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.1 ID 0039140101, KDSS .11533 i 3343 



/i ID onoggosol, ?dss ji336-013O 



7 / 
-// 
















Fig. 3. — For each pair of figures, the left column plots the data for sources with S/N > 4 on a log scale, photons s~^ keV~^ versus 
energy (keV). The data are fit with model D (power-law + intrinsic absorption, data binned to 15 counts per bin). The model, which 
has been folded through the detector response, is plotted as a histogram over the data. Residuals are plotted underneath the fits in units 
of standard deviation. The right column contains T-Nh contour plots, where the Icr, 2it and 3cr contours are plotted as solid lines. The 
plots correspond to sources #:7-#17, starting with #7 in the top left panel. Each pair of plots is labeled with the SDSS name and XMM 
observation ID. 
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Fig. 4. — Oox vs. luv for the red quasar sample, where luv is the 
log of the 2500 A luminosity in units of ergs cm~ s~^ Hz~^. The 
unfilled circles give the measured luv and aox values for the red 
quasar sample. The filled circles give luv,corr and aox,corr, where 
the 2500 A flux has been corrected for E(B- V)spec and the 2 keV 
flux has been corrected for the upper-limit on Njj . The six sources 
with X-ray S/N < 4 are corrected for E(B- V)spec only and not N^f ; 
these sources are marked in red. The S06 olox — luv correlation is 
plotted as a solid line and the S05 sample is overplotted in black 
points. The dotted line marks Ooa: =-1.8. Sources with Cfoa: < —1-8 
are deflned to be X-ray weak. The legend shows the correction to 
aox and luv for an E(B-V)apec = 0.1 and the correction to ciox for 
an N^f upper-limit of lO'^^ cm"'^. Both of the corrections shown in 
the legend were calculated for the samples median redshift, z=1.3. 

= -1.8; those sources with Uox < —1-8 are X-ray weak 
relative to the optical luminosity. All five sources with 2 
< S /N < 4 are X-ray weak, possibly because we were un- 
able to correct the X-ray flux for absorption. All eleven 
sources with S/N > 4 are X-ray normal once corrected 
for absorption (aox.corr > —1-8). 

Correction for absorption can have a large effect on 
F2keV , which can in turn have a large effect on aox ■ The 
largest correction in our sample, for an Njj upper-limit 
of 10^'^ cm~^, led to a factor of 9.5 change in F2keV- An 
additional correction for dust-reddening changed F„„ by 
a factor of 3, leading to A{aox) — 0.2. Sources with 
low S/N and/or high Nh upper- limits will have a large 
but undetermined uncertainty in a^x, which makes it 
difficult to tell if the source is intrinsically weak relative 
to the optical, or if the X-ray weakness is due only to 
absorption. 

There is a weak correlation between aox and T (Spear- 
man Rank = 0.13) such that X-ray weak objects are 
slightly more likely to have flatter spectral slopes. 

4. DISCUSSION 

The redshift selection minimizes the chance that host 
galaxy emission or Lya forest absorption produce the ob- 
served red colors of these quasars. This leaves two pos- 
sible causes of red colors: dust-reddening or intrinsically 
red optical continua. Table [Tl] summarizes the source 
characteristics of the SDSS/XMM red quasar sample, 
which guide the classification of a source as 'absorbed' 
or 'intrinsically red.' We discuss dust-reddened quasars 
in §4.1 and intrinsically red quasars in §4.2. §4.3 dis- 
cusses two quasars that we were not able to classify. We 
comment on the general properties of the red sample in 
terms of the aox-luv relation (§4.5) and the X-ray spec- 
tral slope (§4.4). Finally, we discuss properties of the 
intrinsically red quasars in (§4.5). 

4.1. Dust-reddened Quasars with X-ray Absorption 



By the criteria shown in Table [TTl 8 of 17 quasars in 
the sample are dust-reddened. E{B-Y)spec is significant 
at >3tT for all eight sources. Seven of 8 sources have opti- 
cal continuum shapes consistent with dust-reddened cur- 
vature. (The exception, SDSS J1114+5315 (#11), has 
larger errors on the relative colors, resulting in an am- 
biguous continuum shape. However, strong E(B-V)spec 
and X-ray weakness argue in favor of dust-reddening.) 

SDSS J09 44-H0410 (#1) is a low-ionization BAL quasar 
(|Voit et al.llT993 1. with Si IV, C IV, C III, C II] and 
Mg II absorption troughs of widths ~2,000-5,000 km s^^ 
(Fig.[2|). BALs are kn own to be heavily absorbed in the 
optical and the X-rays (jGreen "eTalllMlh . This source is 
undetected in the X-rays (upper-limit in Fig. [3]). Due to 
the numerous absorption lines and a strong Fell bump, 
we were not able to fit a reddening to the spectral contin- 
uum, so the flatness of the spectrum could be due either 
to line absorption, continuum absorption, or an intrinsi- 
cally different continuum. However, some research sug- 
gests that BALs tend to have redder continua than non- 
BALs, and dust-redde ned templat es appear to be a good 
flt to BAL continua (lYamamoto & Vanscvicius 19991 
Brotherton et al.ll200lUTolea et al...2002. : .Reichard et all 
20031) . We therefore conclude that this quasar is likely 
to have dust-reddening and X-ray absorption. 

All 8 sources match some of the criteria for X-ray ab- 
sorption. For the three sources with 2 < S/N < 4 (where 
X-rays are detected but no X-ray spectral flts can be 
made), all are X-ray weak when the optical luminosity 
is corrected for dust-reddening {aox.corr < —1-8). These 
sources are discussed further in §4.4. 

Four of the 8 dust-reddened quasars have X-ray S/N 
> 4, which allows us to flt a power-law (model C) and a 
power-law plus intrinsic absorption (model D). All four 
sources prefer model C over model D, but two quasars 
have abnormally flat (r<1.5) photon indices to 90% con- 
fldence: SDSS J1226-I-0130 (#8) and SDSS J1435-I-4841 
(#9). This may be a sign of large Nh columns. In 
addition, the upper-limits on N// for all four sources 
are rather high (9 x 10^^ - 1 x 10^3 cm'^, Table l]), 
and are consistent with the gas column predicted from 
E(B-V),pec {NH,opt = 7.0 X 1020 - 1 X 10^1 cm"^) us- 
ing the Galactic gas-to-dust ratio (jBohhn et al.l 119781 : 
iKent et anil99l[ ). Even the ten time s higher Np ex - 
pected from an SMC gas-to-dust ratio (Issa et al."199(^ 
which is more typical for quasars (Maccacaro & Perol^ 
[19811: iMaiolino et al.l[200H I Wilkes et al.ll206^ . is stiU be- 
low the Nh upper-limit found in the X-rays. 

Note that the quasars classifled as dust-reddened oc- 
cupy the lower half of the red quasar sample with respect 
to S/N. The low S/N for these quasars is likely due to 
the photoelectric absorption of the X-rays, as described 
above. The effect of N/f on the X-ray spectrum is further 
discussed in §4.5. 

4.2. Intrinsically Red Quasars 

Seven sources display mild or insigniflcant dust- 
reddening: E(B-V)spec < 0.1 or E(B-V)spec detected 
at < 3cr. Since the reddening in these sources is in- 
distinguishable from a red power-law (§3.2), we classify 
them as intrinsically red using two additional indicators: 
(1) the continuum shape as determined by the relative 
photometry (§3.2) and (2) the lack of X-ray absorption 
(§3.4). 
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TABLE 11 
Summary of Source Characteristics 



ID 


SDSS name 


r 




" Nh,x 


E(B- V)spec 

10^2 em-2 


E(B — y)spec ^ 


Optical 
Continuum 


Notes 


Absorbed Quasars 


1 


SDSS J0944+0410 


1.9 


<-1.69 








Dust 


BAL 


3 


SDSS J1652+3947 


1.9 


-1.92 




0.08±0.02 




Dust 




4 


SDSS J0728+3708 


1.9 


-1.82 




0.25±0.02 




Dust 




6 


SDSS J2217-0821 


1.9 


-2.31 




0.17±0.02 




Dust 




7 


SDSS J1533+3243 


^■^-0.4 


-1.27 


<7.0 


0.12±0.01 


<100 


Dust 




8 


SDSS J1226+0130 


-■■•"-0.4 


-1.66 


<13.4 


O.lSitO.Ol 


<143 


Dust 


NAL 


9 


SDSS J1435+4841 


1 9+0.3 
■■■"^-0.3 


-1.38 


<6.8 


0.13±0.01 


<100 


Dust 


NAL 


11 


SDSS J1114+5315 




-1.66 


<0.9 


0.18±0.03 


<8 


Undefined 




Intrinsically Red Quasars 


5 


SDSS J1002+0203 


1.9 


-1.87 




0.052±0.005 




Red power-law 




12 


SDSS J1133+4900 


1 4+0-5 


-1.52 


<0.8 


0.04±0.01 


<20 


Red power-law 




13 


SDSS J1254+5649 


1 6+01 
^■"-0.1 


-1.30 


<0.2 


10+0" 
"■-^"-0.08 


<3 


Red power-law 




14 


SDSS J0959+0209 


-, 7^+0.09 


-1.46 


<0.05 


0.12±0.03 


<0.6 


Undefined 




15 


SDSS J0958+0213 


-1 70+0. 07 


-1.33 


<0.05 


"■-^"-0.07 


<0.5 


Red power-law 




16 


SDSS J1053+5735 


+0.04 
'^-0.04 


-1.45 


<0.06 


0.11±0.02 


<0.9 


Red power-law 




17 


SDSS J0913+5259 


1 7Q+0.03 
^•'^-0.03 


-1.60 


<0.03 


0-0»0_0.003 


<0.6 


Red power-law 




Unclassified Quasars 


2 


SDSS J0232-0731 


1.9 


-1.95 




0.03±0.02 




Dust 




10 


SDSS J2217-0812 


1 q+0-4 


-1.44 


<1.8 


ii+o-oa 


<20 


Undefined 





Note. — For absorbed quasars, a value is bold-faced if the source matches one of the criteria for dust- reddening in the optical/UV or 
absorption in the X-rays. Criteria for dust-reddening are: a) strong spectroscopic reddening (E(B-V)spec > 0.1) detected to > 3(t confidence 
and b) an optical continuum shape that indicates dust-reddening. The criteria for X-ray absorption are: a) weak X-ray flux relative to the 
optical {aox < —1.8), b) a fiat X-ray spectral slope (F < 1.5 to greater than 90% confidence), and c) a large gas column upper-limit (Nh 
upper-limit > 10^^ cm~^ if intrinsic absorption fits are applied. For intrinsically red quasars, a value is bold-faced if the source matches 
one of the criteria of an intrinsically red power-law: a) reddening is mild {E(B-V)spec < 0.1) or detected at < 3cr, b) an optical continuum 
shape that contra-indicates dust-reddening, c) a low X-ray gas column upper-limit (Nh upper-limits < 10^^ cm~^), and d) gas/dust ratio 
smaller than the Galactic value. 

The optical continuum shape is inconsistent "with dust- 
reddening for all but one object. The relative col- 
ors for SDSS J0959-h0209 (#14) cannot distinguish be- 
tween dust curvature and a red po"wer-law; ho^wever, this 
source meets all of the other 'intrinsically red' crite- 
ria. Six of seven sources also display no significant X- 
ray absorption, with low N^f upper-limits of ~ 10^*^ - 
10^^ cni^^. These same six sources are X-ray normal 
{aox.corr > —1-8), and are not abnormally X-ray flat 
(r>1.5 to 90% confidence). It is not surprising then that 
these six sources have the highest X-ray S/N of the red 
quasar sample, due to the absence of X-ray absorption. 
The exception is SDSS J1002-h0203 (#5), whose low S/N 
(S/N — 2.85) prohibits fitting a power-law or intrinsic ab- 
sorption. This source is X-ray weak {aox ~ -f-82), but 




the low amount of dust-reddening (E(B-Y)spec — 0.052 
±0.005) suggests that the X-ray weakness is intrinsic. 
If we assume an SMC gas/dust ratio, then a reddening 
of 0.05 leads to an X-ray absorption column of Nh = 
3 X 10^^ cm~^. However, a column of less than 10^^ 
cm~ ^ will no t significantly affect the rest-frame 2 keV 
flux llTucker|[l97&) . 

For the four sources with the highest X-ray S/N, SDSS 
J0959-h0209 (#14), SDSS J0958-I-0213 (#15), SDSS 
J1053-h5735 (#16) and SDSS J0913-h5259 (#17), the 
low Nh upper-limits imply gas-to-dust ratio upper-limits 
(Table [8|) up to a facto r of ~2 smaller than the Galac- 
tic value (jBohhn et all 1 19781 iKent et al.l ll99lD and a 
factor of -20 below the SMC value (jlssa et al.l Il990f) 
(Figure [5]). The small dust-to-gas ratio found for these 
four objects contrasts with previous studies, which have 
shown that an SMC-type extinction curve and high gas- 
to-dust ratio is most appropriate for quasar absorp- 



0.4 0.6 
E(B-V),p,, 

Fig. 5.— Nh (lO^o cm-^) vs. E{B - V)spec for the eleven 
quasars with X-ray S/N > 4. Because the Nh column is a 90% 
upper-limit for all of the quasars in the red sample, this results in an 
upper-limit on the gas/dust ratio. Sources classified as intrinsically 
red (§4.2) are marked with a solid circle. The Galactic (solid line) 
and SMC ( dashed line) gas-to-du st ratios are plotted for reference. 
Points from lMaiolino et al.l 1120011 ) are plotted as black points, with 
median error bars shown for one of these points. 

tion (iMaccacaro k Perola I ll98lUMaiohno et al.l l2001l : 
IWilkes et al.ll2002D . For~ample. iMaiolino et'aLr(|200lD 

collected X-ray and optical information for 19 objects 
from the literature. They find gas-to-dust ratios 3-100 
times higher than Galactic for sixteen AGNs. The three 
AGNS with gas-to-dust ratios a factor —1.5 below the 
Galactic value were all low luminosity AGNs {Lx < 10'*^ 
ergs s""'^). For comparison, all of the red quasars have 
> 10^^ ergs s^^, with the exception of the BAL quasar 
(#1). The sma l ler fill ed circles in Fig.[5]mark points from 
IMaiolino et al.| (|2001l). Median err or bars are displayed 
on one of the IMaiolino "eTall (HqoU) points. 
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Fig. 6. — The gas/dust ratio, h ,x-E(B-^) spec, plotted against 
the A(ii — r) — A(g — i) color. The vertical dotted line at A(u — 
r) — A{g — i) = divides the quasars with dust-reddened continuum 
shapes from the quasars with red power-law continua, as described 
in §3.4. The horizontal dashed line shows where the gas/dust ratio 
equals the Galactic value; sources below this line have unusually 
low gas/dust ratios and therefore are considered intrinsically red 
(§4.2). Sources ultimately classified as dust-reddened are plotted as 
open circles, while sources classified as intrinsically red are plotted 
as filled squares. The unclassified source (#10, §4.3) is plotted as 
a star. 

The unusually low apparent gas-to-dust ratios for the 
four high X-ray S/N objects (#14-17) suggest that dust- 
reddening and X-ray absorption are not applicable in 
these cases. We therefore take the unusually low gas-to- 
dust ratios of these four sources to be another indicator 
of an intrinsically red continuum. 

The 'intrinsically red' classification is shown in Fig- 
ure [6l where we plot the gas/dust ratio, Nh,x/E(B- 
V)spec, against the A{u — r) — A(g — i) color. The 
vertical dotted line at A{u — r) — A{g — i) = di- 
vides the quasars with dust-reddened continuum shapes 
(A(m — r) — A{g — i) > 0) from the quasars with red 
power-law continua {A{u — r) — A{g — i) < 0), as dis- 
cussed in §3.4. The horizontal dashed line shows where 
the gas/dust ratio equals the Galactic value. The sources 
in all but the top right quadrant are inconsistent with 
dust-reddening due to their low gas/dust ratios relative 
to the Galactic value and/or the lack of dust curvature 
in their optical continua. 

Table[5]shows that the dust-reddened quasars are more 
likely to have the largest color excesses, while the intrin- 
sically red quasars are more likely to have the smallest. 
This difference is significant in a K-S test at a probabil- 
ity of 1.7% that the two populations come from the same 
parent color distribution. This suggests that, while dust- 
reddening can occur to arbitrarily high values of E(B-V), 
the intrinsic mechanism cannot steepen the continuum 
beyond a critical limit. 

The mean optical slope of the intrinsically red quasar 
candidates is < Uopt > — -1-30, where cc v", with a 
dispersion ctq^pj = 0.34. Figure [7] shows a solid histogram 
of the optical slopes for the intrinsically red quasars, com- 
pared with an open histogram of the optical slopes for the 
Large Bright Quasar Survey (|Francis et al.lll99l[ ). The 
intrinsically red quasars clearly lie in the red tail of the 
LBQS distribution. For comparison, the power-law ob- 
tained for the Vanden Berk composite SDSS quasar spec- 
trum (2001) between 1300 and 5000 A is aopt = -0.44. 

4.3. Unclassified Quasars 



50 r 

40 r 




Fig. 7 . — The optical spec tral indices for the Large Bright Quasar 
survey IIFranc is et al.»199H) are shown in an open histogram. Over- 
plotted is a solid histogram showing the optical spectral indices for 
the intrinsically red quasars. 

The classification of two quasars is undetermined. 
SDSS J0232-0731 (#2) is X-ray weak {aox = -1.95), but 
the signal is too weak to fit F or intrinsic absorption (S/N 
= 3.23). Dust-reddening is mild, and is significant only 
at 1.5(7 (E{B-Y)spec = 0.03 ±0.02). However, the optical 
continuum shape is consistent with dust-type curvature 
(Tabled. 

SDSS J2217-0812 (#10) is also ambiguous because, 
while it has normal F and aox.corr values (F = l.Oj^Q g, 
<Xox.corr = -1.44), the uppcr-limit on Nh is fairly high at 
1.8 X 10^^ cm~^, and the gas-to-dust ratio upper-limit 
is between the Galactic and SMC values. The optical 
continuum shape is consistent with both dust-reddened 
curvature and a red power-law within the photometric 
error bars on the relative colors. While absorption is not 
specifically indicated, it also cannot be ruled out. 

4.4. Red Quasars and the aox-luv relation 

Several studie s fAvni fc Tananbauml 119821: 
Vignali et ai] l2003t IStrateva et all 120051 : ISteffen et all 
20061) ■ have investigated an anticorrelation between aox 
and luv (the log of the 2500 A luminosity) in radio-quiet 
quasars. Figure 2] (described in §3.2) shows that the 
majority of the sample is consistent with the Strateva 
et al. (2005, hereafter S05) Uox - luv sample. Five 
sources, though, are X-ray weak, even after a correction 
for E(B-V)spec is applied, and appear inconsistent with 
the S05 aox - luv sample. All five have S/N<4, so no 
correction can be made for N//. 

We have made a statistical comparison between the 
red quasar sample and the S05 sample. The S05 sam- 
ple was selected from the SDSS using Data Release 2 
photometry, so aside from updates to the SDSS quasar 
selection process between DR2 and DR3, the S05 quasar 
selection is the same as that of the red quasar sample, 
except for the red color cut. We simplified the compar- 
ison of the two samples by noting that the red sample 
luminosity is spread over a single decade, and for S05, 
oiox changes by only 0.16 in this luminosity range. We 
thus restricted the luminosity range of both samples so 
that we could apply the one-dimensional Kolmogorov- 
Smirnov (K-S) test. The observed luminosity range of 
both the red quasar and the S05 samples are restricted 
to luv = 29.4-30.3 when comparing the uncorrected red 
quasar sample (Fig. [4l open circles) to the S05 sample. 
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Since correcting Uox for dust-reddening shifts the UV lu- 
minosity to a brighter interval, we restrict the luminosity 
range to luv = 29.7-30.8 when comparing the corrected 
red quasar sample (Fig. |4l filled circles) to the S05 sam- 
ple. This excludes only one member of the red quasar 
sample (#17). 

First we compare the complete red quasar sample (ex- 
cluding the high Z„„ source #17) with the 805 sample. 
We find a K-S probability of 0.02% that the two sam- 
ples come from the same parent population. Next, we 
consider only the 11 sources for the K-S test where we 
were able to correct for X-ray absorption (i.e. S/N > 
4, black circles in Fig. H]), so that we can compare the 
uncorrected and corrected aox- luv distributions to the 
S05 sample. Comparing first the uncorrected red quasar 
sample to the SOS sample, we find a probability of 10% 
that the two samples come from the same parent popula- 
tion. Then we compare the corrected red quasar sample 
to the SOS sample, and we find a a K-S probability of 
76%. 

Six of 11 quasars with S/N>4 are intrinsically red, so 
this suggests that intrinsically red quasars are not in- 
trinsically X-ray weak compared to typical quasars, with 
the exception of SDSS J1002-h0203 (#5). However, their 
peculiar SEDs make the physical interpretation unclear. 

For sources with S/N<4, it is plausible that correcting 
for absorption would significantly increase the K-S prob- 
ability if the low S/N is due to large Nh columns. We 
find that if we apply a uniform correction to all of the 
low S/N quasars, an absorption column of Nh > 10^^ 
cm~^ is required to achieve a K-S probability of at least 
1% that the red quasar sample is drawn from the same 
parent population as the SOS sample. 

4.5. X-ray Power-law Index 

The mean X-ray photon index for sources with S /N > 
4 is < r > = l.S8±0.09, with an intrinsic dispersion ar 
— 0.3. Figure [8] plots F vs. S/N for each source with 
S/N > 4, with the mean F plotted as a solid line for 
reference, and the intrinsic dispersion plotted as dashed 
lines. For sources with multiple observations, the photon 
indices are averaged and the error plotted is the rnis of 
the measured errors. 

The mean F in this sample is flatter than that of 
other samples of quasar X-ray spectr a including XMM - 
Newton samples. For example, Mat eos et al.l ()200Sf ) 
found a weighted mean < F > = 1.96±0.01 with an 
intrinsic dispersion or = 0.4 for 1137 serendipitously- 
selected XMM-Newton AGN. An investigation of 86 
brigh t AGN d etected in the XMM-Newton GOSMOS 
field (i Mainieri et al. 2003) gives < F > = 2.08±0.08, 
with an intrinsic dispersion ar = 0.24. The two surveys 
are consistent within errors. These mean values are plot- 
ted as dotted lines for reference in Figure [3 

The measurement of Nh can affect the measurement 
of F. The contour plots in Figure [3] show that the two 
variables are correlated so that the same spectrum may 
be fit with a large N// column and steep F, or with a 
smaller Nh column and flatter F. Therefore, the flatter 
mean F may be due to the effect of Nh on the X-ray 
spectra with lower S/N. Because the sources with higher 
X-ray S/N have tighter upper-limits, the F-N/f correla- 
tion is less noticeable for these sources (e.g. Fig. [31 #14, 
15, 16, 17). 
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Fig. 8. — Best-fit X-ray spectral slope (F) vs. signal-to-noise. 
Objects classified as intrinsically red are marked as filled squares, 
while the dust-reddened objects are marked as open circles. The 
unclassified source (#10, §4.3) is plotted as a star. The average F 
and the intrinsic dispersion for the red sample are plotted as solid 
and d ashed lines, respectively. The a verage F for t he Mateos et al] 
BOOSI) sample (< F > = 1.86) and the lMainierTeFaLl (i2007l) sample 
(< F > = 2.08) are shown as dotted lines. 

Gonsidering only the the intrinsically red group 
(squares in Fig. [5]), excluding #5 due its low X-ray S/N, 
we find a mean F = 1.66±0.08, with an intrinsic disper- 
sion of only (Tp — 0.06. This tight intrinsic dispersion 
shows that most of the variance in the intrinsically red 
group comes from measurement errors. Since the contour 
plots in Figure [3] show that F and Njj are not strongly 
correlated for any of the intrinsically red quasars, the 
flatter X-ray spectra cannot be explained by absorption. 

A significantly flatter F may be explained via rela- 
tion to other quantities. The X-ray spectral slope of 
normal quasars does not de pend on optical lurn inosity 
or redshift (jShemmer et al.l 120051: IVignah et ar '2005(1 
but does depend slight ly on X-ray luminosity (^Dai et aT] 
12004 ISaez et al.ll2008D . However, the F-L^ correlation 
contains too much intrinsic scatter (e.g. Figure 8 of Saez 
et al. 2008) to make conclusive statements about the 
intrinsically red quasars. The F-L^, corrrelation is not 
significant in the red quasar sample. 

The strongest correlation reported between X-ray pa- 
rameters and other variables is between F and the H/3 
emission-line FWHM, where flatter X-r ay spectral slopes 
correlate with broad er emission lines (jLaor et al.l 119971 : 
iShemmer et aH [20061) . The H/3 emission line width is in 
turn believed to be anti-correlated with the accreti on rate 
(jBoroson fc Green I [1991 iBrandt fc Boiler l[l998l ). Nar- 
row Line Seyfert I's (NLSl) lie on one end of the F-H/3 
correlatio n, with steep X -ray spectral slopes, narrow H/3 
hues (Bra ndt fc Boiler I 1998: Grupc 2004), and accre- 
tion rates close to the Eddington rate (jKomossa Il2007f l. 
While the H/3 line is not visible in the optical spectra for 
the red quasar sample, the broad Mgll line is visible in 
all of the intrinsically red spectra and has an ionization 
potential similar to H/3, and a FWHM that co rrelates 
weft with FWHM (H/3) (jMcLure fc Jarvis I [2003 ). The 
Mgll line widths of all the red quasars are broader than 
average (Figure [9]) . This may be due to our color and 
redshift selection criteria, which biases the red quasar 
sample towards objects with strong Mgll lines. However, 
the unabsorbed, flat X-ray spectral slopes combined with 
the broad Mgll lines put the intrinsically red quasars on 
the correlation at the opposite end from NLSl's (Fig- 
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equation from iMcLure fc: Jarvis I ()2002f ): 
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Fig. 9.— The histo gram shows the Mgll hne widths (FWHM) of 
the full sample in Sh en et al. I 1)20071 ) (open histogram), the dust- 
reddened red quasars (line histogram), and the intrinsically red 
quasars (filled histogram). 
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Fig. 10.— The wel l-known FWHM(H/3-r relation is shown, us- 
ing the lGrupe I 120041) sample of Seyferts for reference. Broad-line 
Seyfert I's are plotted as open circles and narrow-line Seyfert I's 
are plotted as filled circles. For the red quasar sample, the width of 
the Mgll line is plotted instead of the H/3 due to its availability and 
similarity to H/3 (§4.4). The intrinsically red quasars are plotted 
as filled stars. F is not well-determined for the absorbed quasars 
due to the high upper-limits on N^^, so we do not plot them on 
this plot. 



ure fTU)) ^'^. NLSl's are believed to have accretion rates 
close to the Eddington rate because of their na rrow H/3 
lines and steep X-ray spectra ( Kom ossaH l2007f) . Steep 
red optical/UV continua seem to be an effective means 
of selecting this extreme population. 

4.6. Properties of Intrinsically Red Quasars 

In this section, we investigate the basic properties of 
the intrinsically red quasars: their black hole masses 
(Mbh) and accretion rates (M). These values, and the 
values from which they are derived, are listed in Ta- 
ble [111 The black hole masses are calculated by assum- 
ing virial motion of the broad-line clouds around the cen- 
tral black hole, an assumption justified from reverbration 
mapping (Peterson fc Wandel 1 12000: On ken fc Peterson I 
[2002 ). The FWHM of the broad MgH line can then be 
used to determine the black hole mass using the following 

We do not plot the absorbed quasars on this plot because, 
due to the high N^f upper-limits, F is not well-determined for these 
objects. 
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The luminosity at 3000 A (L3000) gives the radius 
of the BLR as determined by the Rbl^-L re lation ob- 
tained from reverbera tion studies (e.g., Kasp i et al.l 
l2000HBentz et al.ll2006[ ). This relation was derived for a 
normal BBB-dominated SED. Since the BBB appears to 
be missing from the intrinsically red quasars, the ionizing 
continuum may be lower than normal; thi s would result 
in an overestimate of the black hole mass. IShang et al. I 
()2005l ) show that the BBB peaks at --1100 A, so aox 
is a good estimate of the ionizing continuum. Since the 
intrinsically red quasars have aox values within the nor- 
mal range (cf. Steffen et al. 2006), we can assume that 
the Rbl/j-L relation is still valid for the intrinsically red 
quasars. 

The apparent black hole masses for the intrinsically 
red quasars are relatively high, ranging from ~10^ to 
^IQi" Mq (cf. Shen et al. 2007). 

We obtain the accretion rate from the bolometric lumi- 
nosity: Lboi — T]Mc^. We assume an efficienc y 77 = 10%. 
Standard bolometric correction factors (e.g., lElvis et al.l 
[Tool iKaspj et al.l[20Ml are not reliable for the intrin- 
sically red quasars with their atypical SEDs. Instead, 
we integrated the optical power-law {aopt) from 2,500 
to 10,000 A, the nominal UV-to-X-ray power-law [uox) 
from 2,500 A to 2 keV, and the X-ray power-law (F) 
from 2-10 keV to give a rough approximation of the 
bolometric luminosity. Changing the break point be- 
tw een the two power- laws from 2500 A to, e.g. 1150 
A (jShang et al. Il2005[ ). or even to 1000 A, has a small 
effect on h^oi- Since object #17 is a lensed quasar, we re- 
duce its luminosity by the modeled magnification factor, 
M=15 (|Chartaslf20CO . Although we have no mid-IR 
data for these quasars, IR emission accounts for ~40% 
of the bolometric lum inosity for a typical quasar SED 
([Richards et al.l[2006( l. Therefore, if the intrinsically red 
quasars are similar, the IR emission will not change the 
total luminosity, and so the accretion rate, by more than 
a factor of a few. To constrain this estimate, we searched 
the IRAS Point Source Catalog and the Faint Source Cat- 
alog (FSC) for the red quasars, but found no matches. 
Sources in the FSC^** have flux densities greater than 
200 mjy at 60 //m, which gives an upper- limit to AL>,(60 
/im)/Lboi ^ 7.0. This constrains the IR emission to be 
less than an order of magnitude of the bolometric lu- 
minosity. The estimated bolometric luminosities give a 
lower bolometric correction compared to previous work. 
Lfcoj/ALsooo = 2.6 ±1.6 for the intrinsical ly red sample, a 
factor 2.4 lower than that obtained bv , Elvis et al.l ([1994[ ) 
(6.3 ± 3.1). 

Using Lfeoj, we find that the intrinsically red 
quasars have relatively low accretion rates (Fig- 
ure llip . The intrinsically red quasars lie clustered 
between ~0. 001-0. 03 M^dd- In comparison, typical 
accretion rates for optically or UV-selected quasars 
range from M / M^ dd ^ 0-03 - 10.0, with a me- 
dian value of "'^P-6 ([Warner. Hamann fc Dietrichl [20041 : 
[Bonning et al. I [20071 )^ The low accretion rates of the 

http:/ /irsa.ipac.caltcch.cdu/IRASdocs/survcys/fsc.html 
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TABLE 12 

Black hole masses and accretion rates of intrinsically red quasars 



ID 


SDSS name 


FWHM(Mgll) 


Mbh 












km s~-^ 


109 M0 


10''^ ergs 


M0 yr-i 




5 


SDSS J 1002+0203 


10,800 


9.1 


2.5 


0.5 


0.002 


12 


SDSS J1133+4900 


9,100 


3.3 


4.2 


0.7 


0.010 


13 


SDSS J1254+5649 


8,200 


0.68 


1.6 


0.3 


0.019 


14 


SDSS J0959+0209 


23,200 


9.9 


1.6 


0.3 


0.001 


15 


SDSS J0958+0213 


10,500 


1.3 


2.5 


0.5 


0.015 


16 


SDSS J1053+5735 


5,400 


0.76 


3.4 


0.6 


0.035 


17 


SDSS J0913+5259 


6,000 


3.9 


2.6 


0.5 


0.005 
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Fig. 11. — Black hole mass (Mq) is plotted against the accretion 
rate (M© yr~^). The Eddington accretion rate is plotted for a 
reference as a solid line, while correspondingly smaller fractions of 
the Eddington luminosity are plotted as dash-dotted, dashed, and 
dotted lines respectively, as labeled in the plot. 

intrinsically red quasars are in agreement with their po- 
sition on the r-FWHM(Mgll) correlation (Fig.[ini §4.5) 

5. INTRINSICALLY RED QUASARS: PHYSICAL MODELS 

We have found a substantial subset of red quasars 
whose red colors are probably due to intrinsically red 
optical emission in the optical/UV rather than dust- 
reddening. This subset comprises a large fraction (7/17 
~ 40%) of SDSS quasars with {g - r) > 0.5. Note that 
this is based on color selection; selection via relative col- 
ors (e.g. A{g — i) > 0.35) would likely change the fre- 
quency of intrinsically red quasars because dust-reddened 
quasars have a redder color distribution than intrinsi- 
cally red quasars (§4.2). For redshifts z ^ 1 - 1.5, the 
Small Blue Bump reddens the observed g — r colors re- 
sulting in similar samples, whether selected by color or 
color excess. However, for redshifts z 1.5 - 2, the g — r 
color decreases; for this range, selection via color excess 
would increase the number of quasars with smaller color 
excesses, thereby possibly increasing the fraction of in- 
trinsically red quasars. 

Quasars with intrinsically red SEDs have been dis- 
cussed before (e.g., Risaliti et al. 2003; Hall et al. 2006, 
see §1). For e x ample , the intrinsically red quasars found 
by iHall et al.l ()2006l ) do not show absorption in the op- 
tical or X-ray continua, yet they also have steep optical 
slopes (< aopt > — -1.0, Ca^^^ = 0.2), determined via fits 
to the SDSS spectra. 

Intrinsically steep optical power-laws suggest that 
these quasars have a different continuum emission mech- 
anism from ordinary quasars. We consider 3 possibil- 
ities to explain the red optical power- laws: (1) strong 
synchrotron emission visible in the optical, (2) a high- 
temperature BBB, which exposes the red power-law that 



may underlie typical optical emission, or (3) a low- 
temperature BBB resulting from a low accretion rate, 
such that the steep, high-energy tail of the BBB is visi- 
ble in the optical. 

5.1. Synchrotron Emission 

Synchrotron radiation could result in a red optical 
power-law. However, typical synchrotron emission can 
be discounted for two reasons. First, synchrotron radia- 
tion is strongest i n the radio, wit h a turnover frequency 
at = 10^^ Hz (jMarscheilll995| ). Therefore, we would 
expect strong radio emission in all of the intrinsically red 
quasars, when in fact only one can be classified as radio- 
loud (Table [T|). Second, even if a higher than normal 
turnover frequency allowed the synchrotron emission to 
dominate in the optical, the superposition of synchrotron 
emission on a BBB would result in a 'U'-shaped spectrum 
(Fig. 2 of Francis et al. 2000). This is not observed in 
any of the optical spectra. 

5.2. Exposed Underlying Power-law 

A red power-law is someti mes proposed to under- 
lie the BBB in other qu asars (Malka n fc Sargent 1 119821 : 
IWard et al. ]|1987t iLawre nce 2005). If the BBB were 
removed or diminished, this underlying power-law would 
be exposed. One way of "removing" the BBB is to shift 
it to higher temperatures. Possible examples of high- 
temperature BBB's exposing a red, underlying optical 
power-law are kno wn: the narrow-line Seyfert 1 galax- 
ies RE J10 34-H396 (iPuchnarewicz et al.llT995a^ and RE 
J2248-511 (jPuchnarewicz et al.lll995bD . In these objects, 
the optical power-law slopes are -1.3 and -0.9, respec- 
tively RE J1034-H396 is an EUV-selected source with 
an ultrasoft X-ray excess that peaks near 0.4 keV. The 
effective maximum temperature of the accretion disk is T 
^ 10^ K (log Vmax ~16.8), more than an order of magni- 
tude hotter than typical quasar accr etion disk tempera- 
tures (T ~ IQ-*-^ K, log Vmax ~15.3, iMalkan fc Sargent I 
1982). The optical spectrum of this source shows no signs 
of a BBB. RE J2248-511 also has a high-energy turnover 
at around 0.25 keV, but a blackbody is a poor fit to the 
X-ray data. 

The presence of an underlying power-law may be ex- 
plained by weak synchrotron emission from an associ- 
ated jet emitting at a large angle to the line of sight. 
This synchrotron emission would necessarily be weak 
for a radio-quiet object, which is why the BBB must 
be removed or otherwise modified in order for the red 
power-law to be visible in the optical. Alternatively, un- 
saturated Comptonization of a seed spectrum can pro- 
duce a power-law over seven deca des in frequency (e.g., 
iMaraschi. Treves fc Roasio Ill982[ ). 
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Assuming a si mple model o f the accretion disk as a sum 
of blackbodies, iLawrencel ()2005f ) showed that quasar 
SEDs scale homologously with temperature, such that 
hotter AGN have thermal peaks at higher frequencies. If 
the local temperature of the disk is due to the release 
of the gravitational binding energy, then the tempera- 
ture of the accretion disk scales homologously with the 
accretion rate. Since the Mbh determines the inner ra- 
dius of the accretion disk, scaling with black hole mass is 
model-dependent and non-homologous; however, simple 
models show that a less massive black hole will have a 
hotte r maximum effective disk temperature (jLawrence I 
|2005| ). Therefore, high accretion rates and/or low black 
hole masses could shift the BBB to higher tempera- 
tures/energies. 

Neither high accretion rates nor low black hole masses 
are observed in the intrinsically red quasars (§4.6), but 
since the SDSS/XMM red quasars have relatively high 
redshifts (z ~ 1- 2), a soft X-ray tu r nover similar to 
that found in the ' Puchnarewicz et al.l ()1995al |bl) objects 
would be shifted to ~0.4 keV/(l-hz) = 0.1-0.2 keV. This 
turnover would lie below the observed XMM-EPIC en- 
ergy band (0.5-10 keV). A hidden BBB in the EUV or 
soft X-rays would increase the bolometric luminosity, 
thereby increasing the accretion rates. However, an in- 
crease in the ionizing continuum would mean that the 
black hole mass calculations are underestimates. Since 
the effective maximum disk temperature goes roughly as 
T oc (M1/4m-^/2), a NLSl {<Mbh > ~ 10^ - 10^ Mq, 
Grupe & Mathur 2004) can achieve an effective max- 
imum disk temperature of 10^ K at normal accretion 
rates, but highly super- Eddington conditions {h/hEdd ^ 
10^) are required for a quasar with Mbh ^ IO^Mq to 
reach the same temperature. If the black hole masses 
are underestimates, then even higher accretion rates are 
required. Therefore, a high-temperature BBB is unlikely 
to explain the intrinsically red quasars. 

5.3. The Tail of a Low- Temperature Big Blue Bump 

Low disk temperatures due to low accretion rates 
and/or large black hole masses should shift the peak of 
the Big Blue Bump to longer, near-infrared wavelengths 
(Frank, King & Raine 2002, p. 90) . Our bolomet- 
ric luminosity estimates for the intrinsically red quasars 
suggest accretion rates '^10-100 times lower than typi- 
cal quasars (Fig. [Til §4-6). The low accretion rate hy- 
pothesis is supported by the uniformly flat X-ray spec- 
tral slopes and broad Mgll lines of the intrinsically red 
quasars (Fig. [TOl §4.5). We can determine the temper- 
ature of the thermal peak expected for the calculated 
accretion rate u sing the scaling relation developed by 
lLawrenc"e1 ()2005l) . A typical q uasar has an effective max - 
imum temperature -lO"^-^ K (iMalkan fc Sargent I [198^ : 
this corresponds to a BBB peak at log Vmax ~ 15.6. 
Ass uming a co n servat ive accretion rate of O.lMEdd for 
the lElvis et all (fT99l SED, an intrinsically red quasar 
accreting at a rate lower by a factor 10 will be a factor of 
(10)^/'' ^1.8 lower in temperature, so the BBB will peak 
at log I'max ^ 15.3. The most extreme accretion rate 
(0.001 A/^drf, #14) gives a BBB peak at log ly^ax ~ 15.1 

We can model the low disk temperature expected from 
a low accretion rate using a simple template SED. The 
CLOUDY package approximates the AGN SED with a 



toy model (jFerland Il200ll ): 

The A and B constants are normalizations for the opti- 
cal/UV and X-ray terms, respectively. The UV and X- 
ray slopes are power-laws (a„„ and ax), where F^, oc i^". 
The high-energy cut-off, kTc„t, corresponds to the max- 
imum accretion disk temperature. We p erformed a qual- 
itative fit, first to the Elvis et all (| 19941 ) SE P, and then 
to the intrinsically red quasar SEE). For the I Elvis et al.l 
(jl994[ ) SED, we use a„„ = -0.5, = -1, and Tcut = 
10"-^ K. For the intrinsically red SEDs, we change Ux 
to the value obtained in the X-ray fit, and we change 
Tcut by the amount expecte d for the lower acc retion rate, 
calculated according to the Lawrencel ()2005f ) scaling re- 
lation. Since the UV power-law defines the rise of the 
BBB, we do not change a„„. Figure fT2l shows the SED 
data and model fit for all seven intrinsically red quasars. 
The figure demonstrates that a lower disk temperature 
qualitatively reproduces the observed red power-law for 
four objects. Two objects (#15 and #16) are better fit 
with a an accretion rate an order of magnitude smaller 
than that shown. A third object (#5) is not fit well by 
any model. 

6. CONCLUSIONS AND FUTURE WORK 

We have cross-corre l ated t he SDSS DR3 Quasar Cat- 
alog (jSchneider et all l2005f ) with the XMM-Newton 
archive, selecting the reddest {g — r > 0.5), moderate- 
redshift (0.9 < z < 2.1) quasars. We obtain a sample 
of 17 quasars, 16 of which are detected in the X-rays. 
Using both optical and X-ray data to constrain dust- 
reddening and absorption allowed us to distinguish be- 
tween obscured and intrinsically red quasars, although 
two cases remain ambiguous. 

Eight quasars are dust-reddened in the optical and, 
while the X-ray data prefer an unabsorbed power-law, 
the upper-limits are high enough that X-ray absorption 
at the level expected for an SMC dust-to-gas ratio is al- 
lowed. For the three quasars with high enough X-ray S/N 
to fit a power-law -t- intrinsic absorption model, we ob- 
tain upper- limits of 3 - 13 x 10^^ cm~^. For five sources, 
the X-ray S/N is too low to fit spectral models. How- 
ever we can obtain a lower-limit to N/f by comparing 
the an x — liiv distri bution of the red quasar sample to 
that of IStrateva et a l. (2005), a recent study that uses a 
similar selection method, aside from the red color cut. If 
the two samples come from the same parent distribution, 
an absorbing column of at least 10^'^ cm^^ is required for 
the five red quasars. 

Seven quasars display no evidence of X-ray absorp- 
tion, and dust-reddening is contraindicated by the con- 
tinuum shape as determined by the optical photometry. 
These quasars seem to form a group with intrinsically 
red power-laws in the optical/UV. It seems that these 
intrinsically red power-laws may be caused by lower ac- 
cretion rates: low accretion rates are derived from y^BH 
and Lboi estimates, although the M^bh values may be 
biased high by the red SED's. Moreover the unusually 
broad Mgll emission lines place these quasars on the F- 
FWHM(Mgll) relation, also suggesting lower accretion 
rates. 

The intrinsically red quasars are a substantial popula- 
tion with extreme characteristics that are well-suited to 
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Fig. 12. — The SED is plotted as a solid black li ne for the in t rinsic ally red sources #5,#12-#:17, starting with #5 in the top left panel,. 
Units are log(!.'L,y) [ergs s~^] vs. log(i^) [Hz], The lElvis et al.l l |1994l ) SED is overplotted for reference as a dotted black line. A toy SED 
model employing a simple 'disk' shape and parameters as described in §5.3 is overplotted as a solid, blue line to match the Elvis et al] 
II1994I ) SED. A second toy model is calculated, where all the parameters remain the same except that the cut-off temperature is reduced by 
the amount expected from the calculated accretion rate (§5.3). This second model is overplotted as a solid, red line. The optical power-law 
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X-ray follow-up. Since 7 intrinsically red quasar candi- 
dates have been found in the 1% of the SDSS DR3 quasar 
catalog that overlaps with archival XMM-Newton obser- 
vations, as many as 700 intrinsically red candidates may 
exist in the complete catalog. Since the DR6 has more 
than doubled the number of SDSS quasars, we can expect 
1600 intrinsically red quasar candidates, 16 of which have 
been observed screndipitously with the XMM-Newton. 
The intrinsically red quasar candidates in this paper also 
merit further observations in order to confirm the lack of 
dust-reddening and pinpoint the physical cause of the 
red optical power-laws. For example, NIR (e.g. JHK) 
spectra would include Ha and H/3, and so could give an 
independent measure of dust-reddening via the Balmer 
decrement. NIR measurements will also constrain a low- 
temperature disk model, as well as L^oi and so M. With 
optical and X-ray spectra already available, the intrinsi- 
cally red quasar candidates described in this paper are 
well-suited to more complex SED modeling. 
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